Semiconducting Monolayer Materials as a Tunable Platform for Excitonic Solar Cells 
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The recent advent of two-dimensional monolayer materials with tunable optoelectronic properties 
and high carrier mobility offers renewed opportunities for efficient, ultra-thin excitonic solar cells 
alternative to those based on conjugated polymer and small molecule donors. Using first-principles 
density functional theory and many-body calculations, we demonstrate that monolayers of hexag- 
onal BN and graphene (CBN) combined with commonly used acceptors such as PCBM fullerene 
or semiconducting carbon nanotubes can provide excitonic solar cells with tunable absorber gap, 
donor-acceptor interface band alignment, and power conversion efficiency, as well as novel device ar- 
chitectures. For the case of CBN-PCBM devices, we predict the limit of power conversion efficiencies 
to be in the 10 — 20% range depending on the CBN monolayer structure. Our results demonstrate 
the possibility of using monolayer materials in tunable, efficient, polymer-free thin-film solar cells 
in which unexplored exciton and carrier transport regimes are at play. 



of strong electron-phonon coupling, leading to ultrafast 
photoexcited carrier relaxation and consequent thermal- 
ization loss. Despite such inherent material limitations 
and the related constraints they place on the device ar- 
chitecture, polymer and small molecule XSC technologies 
have progressed to impressive power conversion efficien- 
cies, currently up to approximately 11% @,@|. 

Alternative XSC technologies have emerged in recent 
years, including efficient quantum dot based (Io| and 
more recently nanocarbon based XSC [ll|. The key ad- 
vantage of these novel excitonic devices is the possibil- 
ity of altering the HOMO and LUMO levels, the band 
alignment and the optical absorption by using quan- 
tum confinement in nanomaterials rather than differ- 
ent chemistries as in the case of polymers and small 
molecules. Of great relevance is also the possibility of 
novel device architectures [Io| , while retaining the inher- 
ent advantages of solution-based, low temperature man- 
ufacturing on flexible substrates typical of XSC. 

In this context, the recent advent of two-dimensional 
(2D) monolayer materials with tunable optoelectronic 
properties and high carrier mobility — e.g. monolayers 
graphene-BN [TJ, 



Solar cell devices converting energy radiated from the 
sun to electricity have developed into two main families: 
those based on bulk inorganic semiconductors such as Si, 
GaAs, CdTe, and CIGS [lj, in which free charge carrier 
generation follows light absorption without intermediate 
steps, and those based on conjugated polymers and small 
molecules [H, |j| or other materials where following light 
absorption a complex of hole and electron carriers (exci- 
ton) is formed with a binding energy in large excess of 
kT. The latter type, referred to as excitonic solar cell 
(XSC)[|, realizes charge carrier generation by dissoci- 
ating bound excitons at semiconductor heterointerfaces, 
owing to discontinuities across the interface in the elec- 
tron affinity and ionization potential. 

A typical solid state XSC [5| employs a donor-acceptor 
blend of conjugated polymer or small molecule donors 
with high optical absorption in the visible, and fullerene 
derivative acceptors (e.g. C§o, PCBM or similar fullerene 
molecules) @, H, 0]- In such devices, the polymer donor 
dominates (and limits) the key physical processes in the 
solar cell, including optical absorption and transport of 
excitons and charge carriers. In particular, exciton dif- 
fusion lengths of 5 — 10 nm prevent the use of bilayer 
devices, and mobilities lower than 1 cm 2 /Vs limit the 
thickness of bulk hetero junction devices to less than the 
absorption depth (typically 0.1 — 1/xm). 

In addition, tuning the HOMO and LUMO levels of 
conjugated polymers is a challenging task, requiring trial 
and error chemical synthesis of a large number of com- 
pounds; band gaps of less than 1.5 eV are hard to achieve, 
and thus the absorption loss in the red part of the so- 
lar spectrum can be significant Carrier transport 
in conjugated molecules occurs via polarons in a regime 



of graphene [L2|, graphene-BN [ljj |14j, MoS 2 [15Hl7|, 
and other transition-metal oxides and dichalcogenides 
18] — offers renewed opportunities for thin-film XSC. 
Although graphene is a semimetal that can only be used 
in Schottky-type solar cells [l^], other monolayers are 
semiconductors with strong excitonic effects and form ex- 
cellent candidates for XSC. Due to the peculiar nature 
of transport and electron-phonon coupling in monolayer 
materials [20h22| . novel regimes for photovoltaic opera- 
tion can be envisioned, such as hot carrier extraction, 
multiple exciton generation and coherent exciton trans- 
port. Even ultra-thin devices may effectively capture 
sunlight due to the high optical absorption of 2D mono- 
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FIG. 1. (a) CBN monolayer unit cell used in the DFT calcula- 
tions, with an armchair edge between the C and BN domains. 
Shown in the dashed box is a C atom row, with an approx- 
imate width of 0.25 nm. Following the nomenclature used 
here, the unit cell shown is C3(BN)s. Panels (b) to (d) show 
interfaces between a CBN monolayer and (b) PCBM fullerene, 
(c) a (14,0) zig-zag SWCNT, and (d) another CBN layer with 
different composition. 



layers: for example, one layer of graphene absorbs 2.3% 
of the incident intensity in the visible [23j, with a van 
der Waals stacking thickness of only 3.3 A in the direc- 
tion normal to the layer. Tunable optoelectronic prop- 
erties can be achieved by either controlling the structure 
and quantum confinement within the layer, or by stack- 
ing different monolayers to create novel van der Waals 
structures. An example of the former strategy is the 
case of hybridized graphene-BN (CBN) monolayers [l3| , 
whose electronic band gap, optical absorption and exci- 
ton binding energy can be varied by tuning the C domain 
size and shape, due to quantum confinement of excitons 
within the C domains. The mobility in such CBN layers 

synthesized by CVD can be as high as 20 cm /vs mQ. 



Here we show that interfaces between CBN monolayers 
and carbon-based acceptors such as PCBM and semicon- 
ducting single-walled carbon nanotubcs (s-SWCNT) can 
form tunable type-H band alignments [24j . and are thus 
suitable to realize exciton dissociation in 2D monolayer 
based XSC. The power conversion efficiency limit result- 
ing from such CBN interfaces is shown to be tunable, and 
is estimated to be in the 10 — 20% range. Our calculations 
further suggest that even just two stacked monolayers of 
CBN with the proper structure could constitute an ultra- 
thin solar cell with w 3 A thickness. We propose device 
architectures for the experimental study of transport, ex- 
cited state dynamics and power conversion efficiency in 
XSC based on monolayer materials. 



We employ ab initio density functional theory (DFT) 
calculations using the QUANTUM ESPRESSO code j2j| 
on model CBN systems constituted by a monolayer with 
C and BN stripe domains arranged in a 2D superlattice 
and separated by an armchair edge. The CBN repeat 
unit consists of a layer with 8 atom rows and an overall 
composition of C X (BN)( 8 _ X ) , where x is the number of C 
rows in the structure, each 0.25 nm wide (Fig. [TJa)) fiij ]. 

Interfaces containing CBN sheets (Fig. QJb-d)) are 
formed by placing PCBM or s-SWCNT of chirality (10,0), 
(14,0), and (16,0), or a CBN layer at a Van der Waals 
distance of 3.3 A from the CBN monolayer. An or- 
thorhombic simulation cell was adopted, and all struc- 
tures are fully relaxed within DFT to less than 30 meV/A 
in the residual atomic forces. A 15 A vacuum is placed 
in the direction normal to the sheet to avoid spurious 
interactions with the image system. The Perdcw-Burkc- 
Ernzerhof exchange-correlation functional 27) is adopted 
and ultrasoft psuedopotentials [28| are used to describe 
the core electrons. A kinetic energy cutoff of 35 Ry 
was used for the plane-wave basis set and of 200 Ry 
for the charge density, in combination with converged 
Monkhorst-Pack fc-point grids [29[ of up to 24 x 8 x 1. 
For the CBN bilayer calculations, the unit cells consisted 
of two C X (BN)( 8 _ X ) layers with AB stacking in the BN 
domains (B atoms on top of N atoms) [30], and the vdW- 
DF exchange-correlation functional [31| was employed as 
implemented in QUANTUM ESPRESSO. 

For all the interfaces studied here, the DFT derived 
HOMO and LUMO level offsets were estimated as dif- 
ferences in the peaks of the projected density of states 
(PDOS) for the two structures constituting the interface, 
similar to Ref. 32] . Though this method is usually reli- 
able to estimate valence band offsets 33 . accurate con- 
duction band offsets can only be obtained if the DFT 
error on the band gap is comparable for the two ma- 
terials constituting the interface. In order to estimate 
the conduction band offsets with higher accuracy, for the 
CBN-PCBM interfaces we applied GW corrections [33| 
separately to the band gaps of both materials, and the 
GVF+Bethe-Salpeter equation (GW+BSE) formalism to 
compute the optical gap of the CBN monolayers. For the 
CBN layers, both the GW corrections and the GW+BSE 
optical gaps for three compositions (Ci(BN)y, Czt(BN)4 
and C7(BN)i) were taken from our previous work [14j . 
The corrections decrease linearly for increasing C domain 
size, so that corrections at intermediate compositions 
were derived by interpolation. The GW correction for 
the PCBM electronic gap was calculated here using the 
Yambo code 34j within a GqWo update scheme, using a 
plasmon pole model for the dielectric function together 
with cutoffs of 35 Ry and 8 Ry for the exchange and 
correlation part of the self-energy, and up to 1500 empty 
bands. Further details of the many-body calculations are 
reported in the Supporting Material. Fig. Ufa) shows 
the valence and conduction band offsets (AE V and AE C , 
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FIG. 2. (a) DFT valence band offset (AE V , dashed curve) and DFT, GW and BSE conduction band offsets (AE C , solid curves) 
between C X (BN)( 8 _ X ) monolayers and PCBM fullerene, expressed as a function of the C domain size in the CBN layer. AE V 
and AE C are referenced, respectively, to the valence and conduction band edges of the acceptor, as shown schematically above 
the plot. Within our model, AE V is the same at all levels of theory, since the corrections fall entirely on the conduction band. 
The number x of C atom rows in the CBN superlattice structure increases by one unit (from 1 to 7) at each consecutive point 
in the plot, and for each added row the C lateral domain size increases by 0.25 nm. The alignment is consistently found to 
be type-II for C domain sizes of up to approximately 1 — 2 nm depending on the level of theory used, (b) Power conversion 
efficiency contour plot as a function of the CBN donor optical gap and conduction band offset AE C . Constant efficiency level 
curves up to 21% are shown in figure. The diamonds represent the efficiency limits of the three CBN-PCBM combinations 
yielding type-II alignment in the BSE curve in (a). Efficiency values in the 10 — 20% range are predicted depending on the C 
domain size. 



respectively) at different levels of theory for interfaces 
between PCBM and C X (BN)( 8 _ X ) monolayers with dif- 
ferent C domain sizes. Each consecutive point represents 
the addition of one C atom row in the unit cell (Fig. 
[]Ja)), leading to a 0.25 nm increase of the C domain 
size. The DFT band alignment is found to be type-II for 
C domain sizes of up to 1.5 nm, with the PCBM acting 
as the acceptor at the interface, as seen by the positive 
values of AE V and AE C using the convention shown in 
the figure. The trends in AE V and AE C show that the 
HOMO and LUMO energies of the donor can be tuned 
according to the C domain size, yielding unique control 
over the interface band offsets. This results in a tunable 
power conversion efficiency, as explained below. 

To confirm the type-II alignment found within DFT, 
we apply the GW correction separately to the band gaps 
of CBN and PCBM (see Supporting Material). The GW 



corrected AE C values 35], shown in Fig. [Ha) {GW 
curve) yield the same qualitative trends as the DFT re- 
sults, with type-II alignment for C domain sizes of up to 
1.5 nm. 

Assuming that the CBN layer donor is also the main 
absorber in the XSC, the exciton binding energy in the 
CBN layer is the key quantity to determine the possibil- 
ity to transfer a photoexcited electron to the PCBM. To 
address this point, we derived AE C values as differences 
between the LUMO level derived from the optical gap of 
the CBN donor (using GW+BSE) and the GW LUMO 
level of the PCBM acceptor. This combined scheme uti- 
lizing the optical LUMO of the donor and the quasiparti- 
cle LUMO of the acceptor can correctly describe the min- 
imum energy of the exciton formed after photoabsorption 
in the CBN donor, as well as the correct electronic quasi- 
particle level for the addition of the transferred electron 
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to the acceptor. The AE C values derived at this com- 
bined level of theory are shown in Fig. [2ja) (BSE AE C 
curve), and are used below to calculate the power con- 
version efficiencies; within this approximation, the useful 
range for XSC operation is restricted to C domain sizes 
of up to approximately 1 nm, for which the AE C > 
condition is met [361 ]. 

For the three CBN-PCBM cases satisfying this con- 
dition, we estimate a practical limit to the power con- 
version efficiency. The thermodynamic limit to the effi- 
ciency for thermalized carriers and in the absence of non- 
radiative recombination is set by the optical gap of the 
donor through the Schockley-Quisser limit [37|. However, 
efficiency trends and practical limits are far more useful 
for XSC than ultimate thermodynamic limits 38, 39| . 



Next, we analyze the band alignment at CBN- 



Following Scharber et al. 38} , we estimate practical 
values of the maximum power conversion efficiency r\ for 
CBN-PCBM devices with type-II alignment as: 



0.65 • (Ef tA - AE C - 0.3) • J^, P 



J p h(liui) 



d(hu) 



/gopt.d J p h{nuj) d(huj) 

(1) 

where 0.65 is the fill factor (FF), J p h(f^) is the AM1.5 
solar energy flux (expressed in Wm -2 eV _1 ) at the 
photon energy Hui, and E° pt ' d is the optical band gap of 
the CBN donor. 

In Eq. (P), the (E° pt ' d - AE C - 0.3) term is an esti- 
mate of the maximum open circuit voltage (V oc , in eV), 
calculated as the effective interface gap (E° pt ' d — AE C ) 
taken between the HOMO level of the donor and the 
GW LUMO level of the acceptor, minus 0.3 eV, which 
accounts for energy conversion kinetics 38j, |4l| . The in- 
tegral in the numerator is the short circuit current J sc 
calculated using an external quantum efficiency (EQE) 
limit of 100%, while the denominator is the integrated 
AMI. 5 solar energy flux, which amounts to 1000 W/m 2 . 
The efficiency n is thus estimated as the product FF ■ 
V oc ■ J sc normalized by the incident energy flux, in the 
limit of 100% EQE (42 1 . 

Fig. db) shows the efficiency of the three CBN-PCBM 
interfaces with type-II alignment, as a function of the 
CBN donor optical gap and the interface AE C , the lat- 
ter computed as the difference between the CBN optical 
LUMO and the PCBM GW LUMO as explained above. 
For the Ci(BN) 7 , C 2 (BN) 6 and C 3 (BN) 5 monolayers, the 
efficiency values are, respectively, 11%, 15% and 20%. 
The striking efficiency tunability found here is achieved 
by changing the size of the C domain within a 1 nm 
range. Though such small C domains might seem chal- 
lenging to achieve in practice, the immiscibility of C and 
BN in 2D leads to the formation of a large amount of 
sub-nm scale domains, to the point that the single C and 
BN domains cannot be resolved after the monolayer syn- 
thesis [l3| • Strategies for controlling the domain size and 
shape at the atomistic scale in CBN are also being ac- 
tively explored 
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FIG. 3. DFT valence (AE V ) and conduction (AE C ) band off- 
sets at CBN - SWCNT interfaces, shown for combinations 
of the two CBN structures Ci(BN) 7 and C 4 (BN) 4 and the 
three zig-zag nanotubes (10,0), (14,0) and (16,0). AE V and 
AE C are referenced, respectively, to the valence and conduc- 
tion band edges of the acceptor, as shown above the plot, and 
are plotted as a function of nanotube diameter. Both AE V 
(dashed lines) and AE C (solid lines) are positive in all cases, 
indicating type-II alignment. 



SWCNT and CBN bilayer interfaces. Fig. [3] shows the 
DFT valence and conduction band offsets for interfaces 
between the two CBN systems Ci(BN) 7 and C 4 (BN) 4 
(with respectively 0.25 and 1 nm C domain size) and zig- 
zag s-SWCNT with three different diameters. Both AE V 
and AE C are found to be positive (with the convention 
shown in Fig. [3]) for both CBN cases and regardless of 
the nanotube diameter, implying a type-II alignment for 
these interfaces. The band offsets show little variation 
with nanotube diameter, and AE C becomes smaller for 
increasing C domain sizes, similar to the CBN-PCBM 
case. At the interface, the nanotube behaves as the ac- 
ceptor, and should thus be n-doped in a real device for 
optimal performance. Though we do not verify it ex- 
plicitly, the type-II alignment would be retained at the 
GW and GW+BSE levels of theory. For the alignment 
type to be inverted, the GW correction to the nanotube 
band gap would need to be higher than the corresponding 
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FIG. 4. PDOS (left) and HOMO and LUMO orbital iso- 
surfaces (right) for CBN bilayers studied using DFT. The 
energies in the PDOS plot are referenced to the Fermi en- 
ergy Ef. Two monolayer combinations are shown: (a) 
Ci(BN) 7 /C 4 (BN)4 and (b) C4(BN) 4 /C 7 (BN)i. Type-I or 
type-II band alignments are predicted depending on the struc- 
ture and C domain size of the CBN layers composing the 
bilayer. 



correction to the CBN layer band gap, which are, respec- 
tively, 3.25 eV and 1.72 eV for the d(BN) 7 and C 4 (BN) 4 
cases However, only small GW corrections are pre- 
dicted for s-SWCNT in this diameter range (« 10% of the 
band gap (HI), and such corrections certainly cannot be 
as large as 1.7 eV. Similar to the case of the CBN-PCBM 
interfaces, band offset values with a strong dependence 
on the C domain size in CBN are found, which could al- 
low one to tune the solar cell performance by varying the 
structure of the CBN layer. 

The electronic structure of bilayer systems, in 
which the two composing CBN monolayers have different 
atomic structures, is presented in Fig. [4] We show the 
PDOS and the HOMO and LUMO orbitals of two bilayer 
cases among those studied in this work. In both bilayer 
systems, the LUMO localizes within the layer with the 
smaller energy gap. For the Ci(BN)7/C 4 (BN) 4 bilayer 
(Fig. IDJa)), we observe a complete hybridization of the 
valence states, as seen by the perfect overlap of the va- 
lence PDOS, causing the HOMO to delocalize to both 
layers. In contrast, in the C 4 (BN) 4 /C7(BN)i bilayer the 
HOMO is found to be localized on the C 4 (BN) 4 layer, 
i.e. the layer with the larger band gap. In this case, the 



incomplete hybridization of the valence states is seen by 
the slight valence band offset in the PDOS. We found an 
analogous behavior in the Ci(BN)7/C7(BN)i system. 

On this basis, we predict the Ci(BN)7/C 4 (BN) 4 inter- 
face to be a type-I heterojunction with ohmic character 
(due to the absence of barriers for the transport of holes) , 
and the C 4 (BN) 4 /C 7 (BN)i and Ci(BN) 7 /C 7 (BN)i in- 
terfaces to be type-II heterojunctions, in which exciton 
dissociation may be possible within the bilayer. In both 
cases, upon photoexcitation the orbitals are predicted to 
dramatically change their spatial distribution (Fig. 0}, 
which may lead to opportunities for engineering the flow 
of photoabsorbed energy. We note that the alignment 
types found here for the bilayer systems are retained also 
beyond the DFT level of theory, since the GW and the 
GW+BSE corrections to the band gaps are higher for 
the CBN monolayer with the larger band gap 14j. 

Our results thus suggest that in principle it may be 
possible to fabricate a two atomic layers thick (« 3.3 A) 
XSC in the form of a bilayer device, by stacking two CBN 
monolayers with the proper C domain structure. Alter- 
natively, given the tunability of the CBN layer electronic 
structure, a Schottky junction solar cell may be formed 
between a monolayer of graphene and a monolayer of 
CBN. Considering a monolayer absorbance of approxi- 
mately 2% similar to the case of graphene, a structure 
with as few as 30 — 40 monolayers and a thickness of ap- 
proximately 10 nm could suffice to absorb most of the 
incident sunlight with energy above the band gap. 

In addition to the potential technological innovation 




FIG. 5. Proof-of-concept design of a solar cell based on 
2D semiconducting monolayer materials, for the CBN-PCBM 
material combination. Note that the PCBM could be replaced 
by a second 2D monolayer, and the same architecture can be 
extended to other monolayer materials. The shaded oval in- 
dicates the interface where excitons are separated upon illu- 
mination (from the top in figure) , and the arrows indicate the 
direction of carrier diffusion and extraction at the contacts, 
shown as thin metal fingers. 

resulting from the tunability of the parameters of inter- 
est in photovoltaics, XSC based on 2D monolayer materi- 
als could constitute a new platform for the experimental 
study of quantum transport effects (characteristic of 2D 
monolayers) in photovoltaic devices. For example, while 
hot carrier and hot exciton transfer are possible at or 



6 



near an interface [45|, in common XSC they're hindered 
by the diffusion through the bulk of the absorber. How- 
ever, if a solar cell is made with a single monolayer, no 
bulk diffusion is involved and the hot carrier regime could 
be enabled. In addition, the impact of coherent exci- 
ton transport on XSC performance is largely unexplored 
in photovoltaics, due to the fact that exciton transport 
operates in an incoherent regime in reasonably any /xm- 
thick bulk hetero junction solar cell fabricated using stan- 
dard deposition techniques. 

As a test bed for these fundamental effects, we pro- 
pose in Fig. [5] an architecture whereby XSC based on 
2D monolayer materials could be fabricated and charac- 
terized, for example using ultra-fast spectroscopy mea- 
surements [If! to ascertain the presence of hot exciton 
dissociation or hot carrier extraction. 

We remark that while the tunability of the solar cell 
properties presented here originates from changes to the 
composition and domain structure of a single monolayer, 
additional tuning of physical quantities of interest in pho- 
tovoltaics could be achieved by stacking sequences of dif- 
ferent 2D semiconducting monolayers, an alternative ap- 
proach we are currently investigating. 

In summary, we present the idea of XSC based on semi- 
conducting 2D monolayer materials with the potential to 
achieve 10 — 20% power conversion efficiencies in 10 nm 
thick active layers, and show that combinations of CBN 
monolayers and PCBM or s-SWCNT are well-suited for 
the practical implementation of such devices. Even a 
photovoltaic device as thin as two atomic layers of CBN 
(with different C domain sizes, as shown here) or other 
monolayer materials with type-II band alignment holds 
the potential to achieve solar energy conversion at ex- 
ceptionally small length and ultra- fast time scales. The 
unique tunability in 2D monolayer materials of the band 
gap, interface band alignment, exciton binding energy, 
optical absorption, carrier mobility, and electron-phonon 
coupling entails new opportunities for fundamental stud- 
ies and practical implementation of excitonic solar cell 
devices. 
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SUPPORTING INFORMATION 



Details of many-body calculations using the GW and 
BSE methods. 
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